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ABSTRACT 

We report on the structure of the nuclear star cluster in the innermost 0.16 
pc of the Galaxy as measured by the number density profile of late-type giants. 
Using laser guide star adaptive optics in conjunction with the integral field spec- 
trograph, OSIRIS, at the Keck II telescope, we are able to differentiate between 
the older, late-type (~ 1 Gyr) stars, which are presumed to be dynamically re- 
laxed, and the unrelaxed young (~ 6 Myr) population. This distinction is crucial 
for testing models of stellar cusp formation in the vicinity of a black hole, as the 
models assume that the cusp stars are in dynamical equilibrium in the black hole 
potential. In the survey region, we classified 60 stars as early-type (23 newly 
identified) and 74 stars as late-type (61 newly identified). We find that contami- 
nation from young stars is significant, with more than twice as many young stars 
as old stars in our sensitivity range (K' < 15.5) within the central arcsecond. 
Based on the late-type stars alone, the surface stellar number density profile, 
S(/2) oc , is flat, with F = —0.27 ± 0.19. Monte Carlo simulations of the 
possible de-projected volume density profile, n(r) oc r~'^, show that 7 is less than 
1.0 at the 99.73 % confidence level. These results are consistent with the nuclear 
star cluster having no cusp, with a core profile that is significantly flatter than 
predicted by most cusp formation theories, and even allows for the presence of 
a central hole in the stellar distribution. Of the possible dynamical interactions 
that can lead to the depletion of the red giants observable in this survey - stellar 
collisions, mass segregation from stellar remnants, or a recent merger event - 
mass segregation is the only one that can be ruled out as the dominant depletion 
mechanism. The lack of a stellar cusp around a supermassive black hole would 
have important implications for black hole growth models and inferences on the 
presence of a black hole based upon stellar distributions. 
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Introduction 



Over 30 years ago, theoretical work suggested that the steady state distribution of stars 
may be signi ficantly different for clusters with a massive black hole at th e center than those 



without (e.g. iBahcall fc Woljll976l : ICohn fc Kulsrudlll978l : lYoundll980f ). Stars with orbits 



that bring them within the tidal radius of the black hole are destroyed and their energy is 
transferred to the stellar cluster. In the steady state, this energy input must be balanced 
by the contraction of the cluster core, which appears as a steeply rising radial profile in the 
number density of stars toward the cluster center. This radial profile is usually characterized 
by a power law of the form n(r) oc r"''', with a power law slope, 7, that is steep er than that 
of a flat isothermal core. For a single-mass stellar cluster, iBahcall fc WoliI (119761 ) determined 
the dynamically relaxed cusp will have 7 = 7/4. The presence of such a steep core profile, or 
cusp, is important observationally because it may represent a simple test for black holes in 
stellar systems where dynamical mass estimates are difficult, such as in the cores of galaxies. 
The stellar cusp is also important theoretically as it is a probable source of fuel for the growth 
of supermassive black holes and i ts presence is o ften assumed in simulations of stellar clusters 
having a central black hole (e.g.. lMerrittll2005l . and references therein). 



Theoretical work has progressed from the early simulations by IBahcall fc WoliI (119761 ) to 
include many complicating effects, including multiple masses, mass segregation, and stellar 
collisions, on the density profile of stellar clusters in the presence of a supermass ive black 



hole 



.e.g. 



Bahcall fc Wolj Il977l : iMurphv et al.l Il99ll : lAlexander fc HopmanI 120081 ). These 



theories predict cusp slopes within a black hole's gravitational sphere of influence ranging 
from 7/4 > 7 > 7/3 for a multiple mass population to as shallow as 7 = 1/2 for a coUisionally 
dominated cluster core. An important assumption of all cusp formation models is that the 
stellar cluster be dynamically relaxed. Without this assumption, the stellar distribution 
would show traces of the cluster origin in addition to the influence of the black hole. 

The Galactic center is an ideal place to test these theories of cusp formation as it contains 
the nearest example of a supermassive blac k hole (Sgr A*, with a mass M, = 4.1 ± 0.6 x 10^ 
Mq, iGhez et al.ll2008l : iGillessen et al.l 120091 ). Located at a distance of only 8 kpc, the radius 
of the sphere of influence of the Galactic center black hole (~ 1 pc) has an angular scale of 
~ 25" in the plane of the sky, two orders of magnitude larger than any other supermassive 
black hole. In order to plausibly test stellar cusp formation theories, the stellar population 
used to trace the cusp proflle must be older th an the relaxation time, which is of order 1 Gyr 
within the sphere of influence of Sgr A* (e.g., iHopman fc Alexanderll2006l ). In the Galactic 



center, late-type red giants (K to M) are the most promising tracers of the stellar distribution 
since they are > 1 Gyr old, bright in the NIR, and abundant. On parsec size scales, these 
stars dominate the flux, so early seeing limited observations of the surface brightness proflle 
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in the NIR at the Galactic center were successfully used to show th at the structure of the 
nucle ar star cluster at large scales has a density power law of 7 ~ 2 (IBecklin fc Neugebauer 
19681 ). Integrated light spectroscopy of the CO absorption band- head at 2.4 fim (which is 
domii iated by red giants) also confirmed this slope down to within ~ 0.5 pc from the black 
hole (IMcGinn et al.lll989l: iHaller et al.lll996l) . These measurements, as well as integrated 
light spectroscopy from Figer et al. ( 200ol ) of the inner 0".3, found a lack of CO absorption 
in the inner ~ 0.5 pc of the Galaxy. It was unclear at the time whether this represented a 
change in the stellar population with fewer late-type stars in the inner region or a change 
in the stellar density profile. Because the integrated light spectrum is biased toward the 
brightest stars, contamination from a few young Wolf-Rayet (WR) stars (age < 6 Myr), 
such as th e IRS 16 sources identified in the early 1990s located between 1-2" from the black 
hole (e.g., lAUen et al.lll990l : iKrabbe et al.lll99ll ). can significantly impact measurements of 



the underlying stellar density. 

Several approaches have been taken to identify the cause of the differe nce in the CO 



equiva lent width in the integrated light within the inner 0.5 pc of the Galaxy. IScoville et al. 



( I2OO3I ) attempted to remove the effects of the bright young stars on the surface brightness 
distribution through narrow band imaging with NICMOS/HST. They found a surface bright- 
ness that began to flatten and drop slightly in the central 0.1 p c. Since fainter young stars 
can still contribute to the surface brightness, iGenzel et al.l ( 120031 ) used adaptive optics (AO) 
imaging with if (1.6 fim) and Kg (2 /im) filters to measure the stellar density using num- 
ber counts, which they argue more robustly traced the population of red giants. They find 
7 = 2.0 ± 0.1 at i? > 10" and 7 = 1.4 ± 0.1 at i? < 10", which sug gests that there may b e 
a stellar cusp at the Galactic c e nter w ithin the range predicted by iBahcall fc WoliI (119771 ). 
Further work by ISchodel et al.l (120071 ) . which included deeper AO imaging, produced a fit 
of a broken power law to the observed projected surface number density profile of the form 
E(i?) oc R~^, where R is the projected distance from Sgr A* in the plane of the sky, and F 
the projected power law slope. Their best fit broken power law model has a break radius at 
Rbreak = 6".0 ± 1".0 (0.22 pc), with F = 0.75 ± 0.1 outside Rbreak and F = 0.19 ± 0.05 as the 
slope of the stellar cusp insi de Rbreak- This indicat ed that the cusp slope is perhaps flatter 



than the range predicted by iBahcall fc Wolfl (119771 ). However, with only number counts, it 
was still unclear how much bias to the slope was contributed by the dynamically unrelaxed 
young stellar population, which represents a large fraction of the bright stars (K' < 14) that 



have been spectroscopically identified in the inner 1" (e.g.. lGhez et al.ll2003l : iPaumard et al. 
2006h . 



To make further progress, it is necessary to discriminate between the early and late-type 
stars. One recent approach has been through narrow-band imaging to target wavelengths 
where there are measurable differences between spectral features of the two types of stars. 
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For example, the CO bandhead absorption features starting at 2.29 /im are very strong 
for late-type stars (K & M), and absent for early- type stars (O & early B). Using the CO 
narr owband images from the Gemini North Galactic Center Demonstration Science Data 
Set, iGenzel et al.l ( 120031 ) found that the fraction of late-type stars declines between 10" and 
1" from the Galactic center, and thus, may be flatter than expected from number counts. 
However, this narrowband dataset was shallower than that of the broadband images used 
in the number counts measurement s and had troub l e diffe rentiating between the two types 
of stars at > 13. More recently, iBuchholz et al.l (120091 ) took deeper narrowband images 
{K < 15.5) in CO and found an even more shallow slope in the inner region of the cluster. 
Their best fit slope of F = —0.17 ± 0.09 at i? < 6".0 suggests that the late-type stars 
may even be slightly declining toward the cluster center. While number counting with 
narrowband images has greater discriminating power than broadband images, this method 
has limitations. First, this method is strongly dependent on how well the dataset can be 
calibrated based upon existing spectroscopically determined samples, which is hmited to 
K< 13. In addition, limits on the CO equivalent widths depend on how precisely the slope 
of the stellar continuum can be measured with narrow band filters. 

A spectroscopically selected sample would be ideal to separate the late-type, presumably 
dynamically relaxed population from that of the young early-type stars. Until recently, much 
of the spectroscopy taken of stars in the Galactic center has been seeing-limited and limited 
in sensitivity to if < 13 (see lCenzel et al.ll200ol : iBlum et al.ll2003l : iFiger et al.ll2003h . These 
surveys found a drop in the fraction of late-type stars toward the central ~ 10", but because 
they were sensitive only to the brightest giants, the uncertainty remained as to whether 
the fainter stars would show the same trend. Nevertheless, these results combined with 
measurements from imaging all point toward fewer stars in the central 0.5 pc than would be 
expected for a relaxed cusp. 

By combining the advantages of high spatial resolution from AO imaging with the 
discriminating power of spectroscopy, the advent of integral-field spectroscopy (IFU) behind 
AO has potential to provide the most sensitive method of measuring the stellar cusp profile 
in the inner 0.5 pc of the Galaxy. Measurements with the SINFONI IFU at the VLT have 
shown that it is possible to identify young stars at better than 100 mas angular resolution 
( iPaumard et al.l l2006l : iGillessen et al.ll2009l ). Here, we report on a high angular resolution 
spectroscopic survey of the central 0.15 parsec of the Galaxy using laser guide star adaptive 
optics with the OSIRIS IFU at the Keck II telescope. The survey reaches a completeness of 
40% at K' ^ 15.5, two magnitudes fainter than previously reported spectroscopic surveys 
of late-type stars in this region. The observations and data reduction method are described 
in Section [2] and [31 Using these observations, it is possible to separate the presumably 
dynamically relaxed late- type population from the young, unrelaxed population (Section |4]). 
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We find that the surface number density of late-type sta rs is significantly fiatter than the 
range of power laws predicted by iBahcall fc WoliI (119771 ) (Section [5]). This measurement 
rules out all values of 7 > 1.0 at a confidence level of 99.73%. In Section El we discuss 
possible dynamical effects that may lead to the fiat observed slope and implications for 
future measurements of the stellar cusp at the Galactic center. 



2. Observations 



Near-IR spectra of the central 4" of the G alaxy were obtaine d between 2006 and 2008 
using the OSIRIS integral field spectrograph (ILarkin et al.l l2006l ) in co njunction with the 



laser guide star adaptive op tics (LGS AO) system on the Keck II telescope ( IWizinowich et al. 



20061 : Ivan Dam et al.ll2006l ). The laser guide star was propagated at the center of our field 
and for low-order tip-tilt corrections, we used the R = 13.7 mag star, USNO 0600-28577051, 
which is located ~ 19" from Sgr A*. The AO system enabled us to obtain nearly diffraction 
limited spatial resolution (~ 70 mas) in our long-exposure (15 min) 2 /im data sets. All 
observations were taken with the narrow band filter Kn3 (2.121 to 2.229 /im) centered on 
the Br 7 hydrogen line at 2.1661 /im, with a pixel scale of 35 mas and a spectral resolution 
of ~ 3000. With this filter, OSIRIS has a field of view of l."68 x 2."24 per pointing and 
offers sensitivity to both early and late-type stars. Early-type stars from OB supergiants 
to B type main sequence young stars show a prominent Br 7 line, while late-type (K & M) 
evolved giants contain strong Na I lines in this wavelength range. 

With the goal of covering an ~ 8" x 6" area centered on Sgr A*, we divided this region 
into 9 fields (see Figured]). Seven of the 9 fields have now been observed and we will refer to 
these pointings by their approximate orientation with respect to Sgr A* (i.e. GC Central, 
GC East, GC Southeast, GC South, GC West, GC North, and GC Northeast). We obtained 
roughly ten 900 sec exposures (~ 2.5 hr total exposure time) for each field, with the exception 
of the Northeast field, for which we obtained five exposures. For all but the central field, 
the exposures were dithered in a nine position pattern that covered the center, sides and 
corners of a 0".2 x 0".2 square. For the field centered on Sgr A*, the dither positions were 
the four corners of a 1" x 1" square in an effort to keep the 1" x 1" surrounding Sgr A* within 
the field of view of all exposures. The survey region is oriented at a position angle (PA) of 
~ 105 degrees, with the longer direction oriented along majo r axis of the disk of young stars 



rotat ing in the clockwise direction in the plane of the sky (jPaumard et al.ll2006l : iLu et al. 
20091 ). Table [1] summarizes these observations. 
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3. Data Reduction 

Data cubes with two spatial dimensions and one spectral dimension were produced using 
the OSIRIS data reduction pipehne provided by Keck Observatory. The pipeline performs 
dark current subtraction, cosmic ray removal, and wavelength rectification. The cubes were 
not corrected for the effects of chromatic differential atmospheric dispersion because this 
effect is small within the wavelength range of Kn3, shifting the centroid of stars by only 0.1 
spatial pixel between the shortest and longest wavelength channel. 

For each field, stellar positi ons were estirnated by running StarFinder, a point spread 



function (PSF) fitting program (jPiolaiti et al.l 120001 ) . on images made from collapsing the 



data cubes along the spectral dimension. First, a mosaic was made by combining the indi- 
vidual dithers within each field with the OSIRIS pipeline using an average for all overlapping 
spatial pixels. An image was then created from the mosaicked cube by taking the median of 
the spectral channels at each spatial pixel. Next, between 1-3 of the brightest, most isolated 
stars were chosen by hand as PSF reference stars for StarFinder. This PSF was then used 
to determine the centroids of all the stars in the mosaic in an initial pass. After the first 
attempt at identifying sources, we ran an additional iteration of StarFinder to improve the 
estimate of the PSF, which is then used in a second pass to determine centroids. The stars 
identified by this process were cross-checked with the n earest epoch of d eep images taken at 



K' using LGS AO with the NIRC2 imager on Keck II (IGhez et al.ll2008l ) in order to remove 
false detections in the noisier OSIRIS data and to obtain accurate broadband photometric 
measurements for these stars. As a final step, StarFinder was run on each individual cube 
using relative positional priors from the deeper mosaic cubes and the PSF as determined from 
the OSIRIS mosaic to obtain precise OSIRIS positions in preparation for spectral extraction 
from the individual data cubes. 

The spectrum of every star was then extracted from each cube by performing aperture 
photometry for the star at each spectral channel. We used an aperture of radius 2 pixels 
(70 mas) to extract the stellar spectrum and the median of an annulus with inner radius 
at 2 pixels and outer radius at 4 pixels to determine the sky background. For regions with 
very high stellar density, such as within the central arcsecond, the background was manually 
determined by the median of patches of sky free of stars as close to the target star as 
possible. The final spectrum of each star was then produced by averaging the spectra of 
that star from all data cubes from that night. When the signal to noise ratio (SNR) per 
spectral channel between the spectra were significantly different, they were combined using 
an average weighted by the SNR (calculated between 2.212 to 2.218 yum). 

The combined spectra have SNRs ranging from over 100 for bright stars (K' < 13) to 
~ 5 for the faint stars (K' ^ 15.5) that have marginal line detections. Photon noise appears 
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to dominate the SNR of the final spectra and thus our abihty to assign spectral types to 
most of the stars. In a ddition, regions with strong gas emission such as along the mini-spiral 



( iPaumard et al.ll2004l ) and within the central arcsecond can be problematic for identifying 
Br 7 absorption lines from early-type stars. Interstellar Br 7 emission can be quite strong 
and can have several kinematic components, which can lead to a systematic bias in the 
measurement of the Br 7 line in the early type stars if the background gas component is not 
completely corrected. Confusion is also a problem around very bright sources such as the 
IRS 16 stars, where we were unable to extract several K' ^ 13 stars that were identified in 
NIRC2 imagefl 

An AOV (HD 195500) and G2V (HD 193193) star were observed on each night to correct 
for atmospheric absorption line^. The A star, having no intrinsic absorption lines other than 
Br 7 and He I in the Kn3 filter, serves as a good indicator of atmospheric and instrumental 
features. To remove the Br 7 and He I line from the AOV spectrum, we replace the region 
between 2.155 to 2.175 with the same region from the G2V star and divided by the 
empirically observed solar sp ectrum shifted to th e radial velocity of the G2V star (similar 



to the procedure described in lHanson et al.l (119961 )). The A star spectrum was then divided 



by a 10,000 K blackbody to remove the stellar continuum. 



4. Spectral Identification 

Given our survey sensitivity, we expect to be able to detect the early-type stars in the 
form of bright Wolf-Rayet stars, OB supergiants, main sequence young stars down to B2- 
B4V, and late-type stars in the form of K-MIII evolved giants. The most important feature 
of this survey is its ability to differentiate between early and late-type stars. This distinction 
is apparent when comparing template spectra of early type stars with those of late type 
giants, even in the limited Kn3 wavelength range; the young stars have featureless spectra 
except for Br 7 and He I at 2.1661 and 2.1641 fim, respectively, while the late-type stars 
lack these high temperature lines and show a plethora of other atomic and molecular lines, 
most prominently the Na I lines at 2.2062 and 2.2090 fim. 

In order to assign spectral types to the stars, we separate all the extracted spectra by 
visual inspection into five groups: 1. Spectra with prominent Br 7 lines and no Na I lines 



-'^NIRC2 is less sensitive to confusion because the plate scale is 10 mas compared to 35 mas used in the 
OSIRIS observations and the higher strehl ratios obtained from the much shorter exposures (30 s vs. 900 s) 

^These stars were observed using natural guide star mode with 30 s integration time and 5 frames per 
star. 
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are identified as early-type stars (e.g., SO-2 in Figure [2]). 2. Spectra with prominent Na I 
lines, but no Br 7 line are identified as late-type stars (e.g., SO-13 in Figure |2]). 3. Low 
SNR (< 5) spectra with no apparent lines are classified as unknown. 4. Bright stars with 
high SNR (K' < 14, SNR > 40) spectra with no late-type or Br 7 lines are identified as 
young. These are most likely late 0-type stars, which have very sma ll or undetectable Br 7 



absorption lines at the resolution of this survey (IHanson et al.l 119961 ) (e.g., S2-74 in Figure 
[2]). 5. Spectra with both Br 7 and Na I lines are identified as late- type if there is strong gas 
emission nearby, which would explain the Br 7 line as an over- subtraction of the background 
gas (2 stars showed these features). In the survey region, we classified 60 stars as early-type 
(23 newly identified) and 74 stars as late-type (61 newly identified). The late and early-type 
stars are reported in Tables [3] and H] and labeled in Figure [H The stars with successfully 
extracted spectra, but have undetermined spectral types are listed in Table O Our ability 
to classify the stars is strongly dependent on the SNR of the spectra, with 100% of the stars 
having SNR > 40 being identified (see Figure [3]). 

Equivalent widths for the Na I lines in the spectra of the late- type stars were measured by 



norma lizing the spectra and integrating between 2.2053 to 2.2101 /im as in lForster Schreiber 



(120001 ). after shifting the spectra to rest wavelengths. The Na I equivalent widths range from 
1.2 to 4.9 A, consistent with that of Kill to Mill evolved giants, with RMS uncertainties be- 
tween 0.1 to 2 A. There does not appear to be a significant correlation between the equivalent 
width and K' magnitude, suggesting that the Na I lines are not able to further differentiate 
between K and M spectral types at the SNR of these measurements (see Figure H]). The 
equivalent widths of the Br 7 lines are measured by fitting a Gaussian to the line and inte- 
grating the model Gaussian parameters. While the equivalent widths of the Na I lines do not 
change with brightness, the young stars show a significant decrease in the equivalent width 
of Br 7 with increasing brightness as shown in Figure HI These measurements support the 
assumption that the bright stars {K' < 13) with no spectral features in the Kn3 wavelength 
range are young stars. 

This survey is complete to about K' = 12 and falls off to about 40% completeness at 
K' = 15.5, measured relat ive to the stars d etected in the same region from imaging, which 



is complete to K' = 16.0 (iGhez et al.ll2008l ). We were able to determine spectral types for 
all stars with extractable spectra and K' < 14.0. The sample is not complete to K' = 14.0 
because the spectra of some bright stars {K' ^ 13) detected in broadband images were not 
extractable due to confusion with the IRS 16 sources. The K' luminosity function of all the 
stars in the survey is shown in Figure El with a comparison to that derived from imaging 
at K'. The K' luminosity functions of the stars in the individual fields are similar to that 
of the total sample, with the notable exception of the western field, where there are fewer 
stars detected both in imaging and in spectroscopy. This difference in number of detected 
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stars is most likely due to poor seeing on the night of observation and higher extinction in 
that region. Figures [7] and [8] shows the locations of the stars with their spectral types on a 
collapsed image of each of the OSIRIS fields. 



5. Results 

With spectroscopic identification of the stars brighter than K' ~ 15.5, we are able to 
separate the presumably dynamically relaxed old stars from the unrelaxed young population. 
The stellar number density profile is constructed in radial bins of 0."5 over the area that 
was sampled in this survey out to a radius of ~ 3", with error bars scaling as a/TV, where N 
is the number of stars in each bin. It is clear that, while the number density of early-type 
stars increases quite rapidly toward the central arcsecond, the late-type stars have a very 
flat distribution (Figure [9]). Outside the central arcsecond, the projected number density of 
early-type stars drops to about half that of the late-type stars. For a quantitative comparison 
to stellar cusp models, we fit a power law to the surface number density of both populations. 
The late-type population is best fit with a power law slope Void = —0.12 ± 0.16, while for 
the early-type stars, Tyoung = 1-51 ± 0.21. For comparison, we also measured the projected 
stellar number density over the same survey region using LGS AO K' imaging from July 2008, 
which is much deeper {K' < 20), but with no spe ctroscopic discrimina tion. We find that 



^imaging = 0.19 ±0.06, Comparable to that found by lSchodel et al.l (120071 ) using only number 
counts. The spectroscopically determined late-type slope is flatter than that determined 
from number counts alone, but is consistent within the ~ 2a formal fitting errors. 

The survey has complete azimuthal coverage out to a radial distance of ~ l."5 from 
Sgr A*. Beyond this radius, the radial coverage extends out to ~ 4", but the azimuthal 
coverage is complete due to two fields missing in the northwest and southwest. We find that 
these missing two fields do not significantly impact the radial profile of the late-type stars, as 
removing one or two of the outer observed fields from the analysis leads to the same power 
law slope, within the measurement uncertainties. In addition, previous obser vations of the 



late-type stars in this re gion have found the distribution to be isotropic (e.g. iGenzel et al. 



20031 : iTrippe et al.ll2008[ ). 



Potential systematic errors to the radial profile measurement may arise from variable 
extinction and incompleteness between each field. In order to study how these two factors 
impact our result, the radial profile was computed using stars with extinction corrected 
magnitudes only over regions that are at least 30% complete. To do so, the extinction 



map reported by iBuchholz et al.l (120091 ) was used to correct for the variations in extinction 



between each star by adding — 3.0 to each star to bring them to the same canonical 
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Ak = 3.0 extinction to the Galactic center (see Table [2] for the average and RMS value 
of Ak in each field). Then, the K' luminosity function was recomputed for each field to 
determine the K' magnitude bin beyond which the completeness falls to less than 30%. 
Table [2] shows this magnitude bin for each field before and after extinction correction. This 
separates the fields into three groups in terms of their 30% completeness : all fields are at 
least 30% complete down to K' = 14.5; all fields except W are complete to K' = 15.5; the 
central field is complete to K' = 16.0. We calculate the radial profiles in each magnitude 
bin using only fields that are at least 30% complete at that magnitude. For all stars with 
K' < 14.5, we use all the fields, while for stars with 14.5 < K' < 15.5, the western field 
was dropped from the measurement; the 15.5 < K' < 16.0 magnitude bin was not included 
because only the central field is complete to at least 30% at that magnitude. The radial 
profiles from the K' < 14.5 and 14.5 < K' < 15.5 magnitude bins were then summed to 
produce the final measurement of the surface density, with errors added in quadrature (see 
Figure [To]). A total of 32 early-type and 63 late-type stars are included in the final radial 
profile measurement. The early-type stars have Tyoung = 1-22 ± 0.16, while for the late-type 
stars. Told = —0.27 ± 0.19. Neither of these values differ significantly from the case without 
extinction and completeness correction. The rest of the analysis will use these values for the 
measured surface number density profile. 

Beyond a radius of 1", there are more young stars detected in the eastern and south- 
eastern fields compared to the other observed regions. This excess is in the direction of the 
major axis of the proje cted coherent d isk of young stars orbiting in the clockwise direction 



in the plane of the sky (iLu et al.ll2009l ). however, after accounting for the completeness and 
differences in Ak, neither the number of young, nor old stars shows a statistically significant 
excess. 



6. Discussion 

The effects of projection of a three-dimensional cluster onto the plane of the sky must 
be taken into account for quantitative interpretations of the measured surface density profile. 
The projection onto the sky plane of a spherically symmetric cluster can be done using the 
integral: 

°° rn{r)dr 



S(i?) =2 / (I] 



R 



where R is the distance from the center of the cluster in the plane of the sky, r is the 
physical radius of the star cluster, S(_R) is the projected surface density, and n(r) is the 
volume number density. Because this integral diverges for all 7 < 1.0, where n{r) oc r~'^, 
the density in the outer region of the cluster must fall off steeply to accommodate a shallow 
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power law at the center of the cluster. While we do not observe a break in the cluster 



density profile within our 4" survey, one was observed by ISchodel et al.l (120071 ) and others at 
scales larger than our survey. In order to account for the outer region of the star cluster, we 
therefore model the number density profile as a broken power law: 



n{r) oc 



71 7" < f break) 



(2) 



Using the measurements from ISchodel et al.l (120071 ). we set rhreak = 8" and 72 = 1.8 (the 
deprojected values of Rbreak = 6" and V = 0.8, from the surface density measurements using 
this broken power law model) to constrain the parameters of the outer cluster density profile. 
We assumed that contamination from young stars is less severe in the region outside of our 
survey, so using an outer cluster density profile measured from total number counts should 
introduce only a small bias. For example, the densit y of bright y oung stars falls off rela tively 
steeply as R''^ ( IPaumard et^boOfil : ILu et al.lboOfll . but see also ISuchholz et al.l koO^h ). To 
determine the constraints on the physical density profile in the inner region, we performed 
Monte Carlo simulations of clusters with 71 between -2.0 and 2.0. For each value of 71, 
10^ realizations were performed. In each realization, 74 stars (the number of late-type stars 
used to measure the radial profile) were drawn from the broken power law distribution after 
correcting for the limited area observed in this survey. The locations of these stars were then 
binned and fitted the same way as the data. We find that for —2.0 < 71 < 0.5, the projected 
inner radial profiles are flat with F ^ (see Figure [TTj) . Figure [12] shows the relationship 
between 71 and F as determined from the MC simulations. Using this relationship and the 
observed inner radial profile, F = —0.27 ±0.16, we can set an upper limit of 1.0 to the value 
of 7i the 99.73% confidence level. 

Although the slope measurement in this survey cannot constrain whether there is a 
'hole' in the distribution of late-type stars or just a very shallow power law within the 
central 4", th e inferred slope is signi ficantly flatter than the range of 71 between 7/4 and 3/2 
predicted by iBahcall &: Wolj (119771 ). In fact, a flat core density proflle without a cusp can 
flt the data equally well. Below, we discuss several plausible dynamical scenarios that may 
deplete the number density of the late-type giants observed in this study, as well as possible 
observational prospects for making further progress. 



6.1. Mass segregation 



The true population of dark mass that consists of stellar remnants such as neutron 
stars or stellar black holes is unknown, but theoretically, there is a strong case for a large 
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population of such remnants due to mass segregation (IMorrisi Il993l ). As these dark rem- 
nants migrate inward, they will scatter the lighter stars outward. Recent simulations by 
Alexander fc Hopmaru (120081 ) showed that a population of ~ 10 Mq stellar remnants can 
sink into a much steeper density distribution (2 < 7 < 11/4) than the lighter stars (3/2 < 
7 < 7/2), if the stellar remnants are relatively rare compared to stars. This results in a very 
dense cluster of stellar remnants, but the predi ction of the stellar cus p from strong mass 
segregation does not differ substantially from the iBahcall fc Wolj (119771 ) values. Thus, mass 
segregation cannot be the only mechanism responsible for the inferred lack of a cusp in the 
evolved red giant population. 



6.2. Envelope destruction by stellar collisions 



The M and K type giants in this survey have radii that vary between 20 to 200 R©, 
which results in s ubstantial cross sections for collisions wi th other stars in this high density 
environment (e.g. iBailey fc Davieslll999l : iDale et al.ll2009l ). This may lead to a preferential 
depletion of red giants toward the center, which would result in a biased measurement 
of the underlying stellar cusp slope, since most of the surviving stars in the cusp would 
be unobserved main sequence stars with smaller radii. The depletion of giants has been 
suggested in the past based upon low spatial resolution spectroscopy of the CO band-head 
at 2.29 jUm jMcGinn et ah 1989 ) as well as in narrow band imaging ( Genzel et al. 2003 



Buchholz et al 



2009h . Theoretically, a coUisionally dominated cusp can be as shallow as 



7 = 0.5, which is within the range of slopes allowed by our upper limits, for an isotropic 
distribution, (i.e., for a distribution function that depends only on e nergy), although thi s 
result no longer holds in the case of a more general distribution function ( iMurphy et al.lll99ll ). 
While the collisional destruction of giants can operate efficiently within the region observed 
in this survey, the effectiven ess of this raechan ism drops very quickly with distance from 
the black hole. For example, iDale et al.l (120091 ) found that collisional destruction of giants 
becomes less likely beyond about 0.1 pc from Sgr A*. While this survey cov ers a region out 
to ~ 0.16 pc in projection from Sgr A*, previous observations such as from iBuchholz et al. 



(I2OO9I ) suggest that the flattening in the late-type giant radial profile extends out even 
further to ~ 0.24 pc in projection. If the assumptions about the steep distributio n s of st ellar 
mass black holes and other stellar remnants in the simulations of iDale et al.l (120091 ) are 
correct, then this would suggest that the collisional destruction of giants is not the dominant 
mechanism for clearing the cusp of stars. Observational constraints on the effectiveness of 
this mechanism can be made by comparing the radial density profiles of giants of varying 
stellar radii. Collisional destruction should preferentially remove stars having larger radii at 
a given distance from the black hole (see Section 16.41 for further discussion) . 
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6.3. IMBH or binary black hole 



A more dramatic change t o the stellar cusp can be achieved through the infall of an 
IMBH. iBaumgardt et al.l (120061 ) performed N body simulations of the effects on the central 
star cluster with an initial cusp slope of 7 = 7/4 due to the infall of IMBHs of mass 3 x 10^ 
M© and lO'^ Mq. They find that the IMBH will deplete the cusp of stars, which will cause the 
inner portion of the cluster to resemble a core profile after the IMBH has spiraled into the 
central black hole, and that it would take over ~ 100 Myr after the IMBH infall to replenish 
the stellar cusp. For the case of a 10'^ M© IMBH, the density profile has a fiat core, with a 
core radius of ~ . 2 pc, which is roughly consistent with that observed in this survey and 



bv lBuchholz et al.l (120091 ). 



The existence of a binary massive black hole in the Galactic center can also contribute 
to the loss of stars in this region. A star passing within the orbit of the companion massive 
black hole will undergo dynam ical interactions with the bin ary which can result in the 
star achieving ejection velocity (IMerritt fc Milosavljevial2005l ). In addition, similar to that 
of an infalling IMBH, a merger with a massive black hole system more recently than the 
dynamical relaxatio n time will result in a system that is still out of equilibrium at this time 
(|Merritt et al.ll2007h . 



Current measurement s impose only modes t limits on the possibility of an infalling IMBH 
and a binary black hole. iTrippe et al.l (|2008l ) found no obvious kinematic signatures of a 
recent infall of an IMBH using 3D velocity measurements of the late-type stars, however, a 
more complete study of the phase space distribution of the stars will be necessary to asc ertain 
the likelihood of an IMBH clearing the cusp of stars (see iGualandris &: MerrittI (|2009l ) for a 
more extensive discussion of current constraints on the existence of an IMBH at the Galactic 
center). Current measurements of the reflex motion of the black hole at radio wavelengths 
place the limit to th e mass of a companion bla ck hole at < 10^ M© with semi-major axes 10'^ 
AU < a < 10^ AU (IReid fc Brunthaleiil2004l ). Precise measurements of the possible reflex 
motion of the black hole using stellar orbit s can provide str onger constraints in the future 
on the presence of a companion black hole (iGhez et al.ll2008l ). 



6.4. Prospects for future observations 

Given the current state of measurements of the cusp density, it is unclear which, if any, 
of these dynamical interactions are responsible for the shallow inferred cluster slope. The 
strong degeneracy between the surface number density and the de-projected volume density 
makes it difficult to determine the extent of the depletion of red giants within the inner 0.1 
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pc. Limits to the line of sight distance of these stars from Sgr A* will be necessary to resolve 
this degeneracy. Within the central arcsecond, it may be possible with direct orbit fits to 
measure the orbital elements of stars located physically near the black hole. For example, at 
least one of these late-type stars, SO-17 (which passed within 65 mas of Sgr A* in the plane 
of the sky in 2007) has an orbit measu rement that places it s periastron distance closer than 



8 mpc (1600 AU) from the black hole (IGillessen et al.ll2009l ). This indicates that while there 



are fewer late-type stars expected in the central region, it is not entirely devoid of these stars. 
Even for stars without well-determined orbits, acceleration limits can be used to constrain 
the 3D distribution of the late-type population, and thus help determine whether there is a 
'hole' or simply a flattening of the radial density profile toward the center of the cluster. 

Determining which mechanism is responsible for the shallow slope of the late-type stars 
is important because it is related to whether we can use these stars as tracers of the under- 
lying stellar cluster. The main sequence, > 1 Gyr old stars, represents a significant fraction 
of both the mass and number density of the stellar cusp, but is currently unobservable spec- 
troscopically; spectral types later than A have K' > 19, compared to our current sensitivity 
limit of 15.5 at K' . Fortunately, in a relaxed population, the red giants should have the same 
phase space distribution as main sequence stars of comparable mass, so their radial density 
profile should be identical. However, the red giants have substantially larger radii than the 
main sequence stars, so they are more likely to be destroyed if collisions are important. It 
may be possible to observationally constrain whether collisions are the dominant mechanism 
responsible for the observed slope of the radial density profile by spectral typing the late-type 
stars more precisely than was possible in this survey. For example, the equivalent width of 
the CO band head at 2.29 /im, which varies over a wider range than the equivalent width 
of Na I observed in this study, should permit one to distinguish between individual spectral 
sub-types. With a factor of 6 ch ange in the radii of red giants between spectral types K3III 



and M5III (Ivan Belle et al.lll999l ). a difference in density profile between different spectral 



types would provide a test of whether collisions affect the distribution of these stars. 



7. Conclusions 

We completed a survey using high angular resolution integral field spectroscopy of the 
inner ~ 0.16 pc of the Galaxy in order to measure the radial profile of the late- type stars in 
this region. The survey reached a depth of K' < 15.5 and is able to differentiate between 
early and late-type stars. The survey provided spectral types for 60 early-type and 74 late- 
type stars, with 23 and 61 previously unreported early and late-type stars, respectively. We 
find a late-type stellar surface density power law exponent F = — 0.27 ± 0.19, which limits 
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the volume number density profile slope 7 to be less than 1.0 at the 99.73% confidence level, 
and even allows for the presence of a central drop in the density of late-type giants. We 
present three scenarios that may lead to the depletion of late-type stars, but are unable to 
constrain the candidate mechanism sufficiently to determine which is dominant. Given the 
current measurements, we cannot yet determine whether the distribution of observed giants 
is representative of the distribution of stellar mass. Being able to infer the underlying dy- 
namically relaxed stellar population will be crucial in order to establish whether the Galactic 
center is lacking the type of stellar cusp long predicted by theory. Obtaining an unbiased 
measurement of the stellar distribution is important because such cusps have considerable 
impact on the growth of massive black holes as well as on the evolution of nuclear star clus- 
ters. Progress in achieving this goal will be possible with improved kinematics and spectral 
coverage in order to break the degeneracy in the surface number density profile to better 
establish the three-dimensional distribution of the stellar cluster. 
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RA Offset (arcsec) 

Fig. 1. — The currently surveyed region overlaid on a K' (2.2 /im) image of the Galactic 
center taken in 2007. Sgr A* is marked at the center with a *. Spectroscopically identified 
early (blue) and late-type (red) stars are marked with circles. Each field is enclosed by 
dotted lines. Some spectral identifications are outside of the marked lines because they were 
found at the edge of the dithers. 
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Fig. 2. — Example of observed spectra in the Kn3 narrowband wavelength region of early 
and late-type stars. The Br 7 and Na I lines are used to differentiate between early and 
late- type stars, respectively. SO-2 is a, K' — 14.1 early- type star with a strong Br 7 line. 
S2-74 is a i^T' = 13.3 early-type star showing a featureless spectrum in the OSIRIS Kn3 filter 
wavelength range. SO-13 is a A'' = 13.5 late-type star showing the Na I doublet. See Figures 
2a, 2b, and 2c in the online supplements for spectra of all stars with spectral types identified 
in this survey. 
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Fig. 3. — The relationship between the SNR and K' for stars observed in this study and their 
spectral type. The spectral types of every star with SNR > 40 was determined successfully. 
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Fig. 4. — Top: equivalent width of the combination of Na I hnes at 2.2062 and 2.2090 /im 
as a function of K' m a gnitu de. Also plotted are equivalent width measurements of Na I 
from iForster Schreiberl (120001 ) with spectral types corrected for the distance and extinction 
(Afc=3.0) toward the Galactic center (red). Bottom: e quivalent width of th e Br 7 hydrogen 
line at 2.1661 /im. Equivalent width measurements from lHanson et al.l (119961 ) are also plotted 
for O and B type stars, corrected for the distance and extinction (Afc=3.0) to the Galactic 
center (red). 
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Fig. 5. — The K' luminosity function from the spectroscopic survey in comparison to that 
found from imaging. In the top plot, we also include for reference, the rough spectral types 
expected to be observable at each luminosity bin assuming Ak = 3.0 and a distance of 8 
kpc. The axis on the right shows the H-K color associated with each of the spectral types 
( iDucati et al.ll200ll : Martins fc Plezll2006l ). The WR stars occupy a range in both K' and 
H-K colors. 
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Fig. 6. — The K' luminosity function from the spectroscopic survey in comparison to imaging 
in each individual pointing. 
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Fig. 7. — Images from collapsing the OSIRIS data cubes along the spectral dimension for 
each individual pointing of the survey. The images are oriented with north up and east to 
the left. Spectroscopically identified early (blue) and late- type (red) stars are marked with 
circles. 
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Fig. 8. — Images from collapsing the OSIRIS data cubes along the spectral dimension for 
each individual pointing of the survey. The images are oriented with north up and east to 
the left. Spectroscopically identified early (blue) and late- type (red) stars are marked with 
circles. 
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Fig. 9. — Plot of the surface number density as a function of projected distance from Sgr 
A* in the plane of the sky for different populations: old (late-type, red), young (early- type, 
blue), and total number counts from K' imaging. 
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Fig. 10. — Plot of the surface number density as a function of projected distance from Sgr 
A* in the plane of the sky for different populations: old (late-type, red), young (early- type, 
blue), and total number counts from K' imaging. These radial profiles have been corrected 
for completeness and extinction using the method detailed in Section [51 
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Fig. 11. — Left: broken power law density profiles with break radius, rhreak — 8-0" and outer 
power law 72 = 2.0, and varying inner power laws 71. Right: the projected surface number 
density profile of each of the broken power laws. The fitted inner surface density power law 
r is fiat for 71 < 0.5. 
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Fig. 12. — Left: the results of the Monte Carlo simulations of different inner radial power 
laws 7 and the resulting measured projected power law of the radial profile, F = — 0.27±0.19. 
Center: using our measurements, we can constrain some of the parameter space for F vs. 
7. Right: we marginalize over F to determine that our measurement constrains 7 to be less 
than 1.0 at a 99.73% confidence level. 
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Table 1. Summary of OSIRIS observations 



Field Name 


Field Center^ 


Date 


Nf ^ 


FHWM= 




(") 


(UT) 




(mas) 


GC Central (C) 


0, 


2008 May 16 


11 


84 X 85 


GC East (E) 


2.88, -0.67 


2007 July 18 


10 


85 X 70 


GC South (S) 


-0.69, -2.00 


2007 July 19 


10 


73 X 63 


GC West (W) 


-2.70, 0.74 


2007 July 20 


11 


110 X 86 


GC Southeast (SE) 


1.67, -2.23 


2008 June 03 


11 


68 X 63 


GC North (N) 


0.33, 2.01 


2008 June 07 


7 


102 X 85 






2008 June 10 


5 


75 X 70 


GC Northeast (NE) 


2.55, 1.27 


2008 June 10 


5 


74 X 68 



^RA and DEC offset from Sgr A* (RA offset is positive to the east). 
''Frames have an integration time of 900 s. 

'^Average FWHM of a relatively isolated star for the night, found from a 
two-dimensional Gaussian fit to the source. 



Table 2. Results and survey completeness 



Field 






At'' 


RMSa;, 


Observed 




Ak 


Corrected'^ 




Completeness 


K 


1 


Completeness 


C 


100 


28 


3.24 


0.04 


15.5 - 16.0 


0.4 


15.5 - 


16.0 


0.4 


E 


37 


25 


3.05 


0.06 


15.5 - 16.0 


0.3 


15.0 - 


15.5 


0.3 


SE 


60 


25 


3.29 


0.10 


15.5 - 16.0 


0.4 


15.0 - 


15.5 


0.5 


S 


35 


6 


3.28 


0.05 


15.5 - 16.0 


0.6 


15.0 - 


15.5 


0.6 


w 


35 


19 


3.44 


0.12 


14.5 - 15.0 


0.7 


14.0 - 


14.5 


0.4 


N 


40 


14 


3.39 


0.12 


15.5 - 16.0 


0.4 


15.0 - 


15.5 


0.4 


NE 


32 


7 


3.13 


0.07 


15.0 - 15.5 


0.5 


15.0 - 


15.5 


0.5 



^Number of stars detected by StarFinder. 
''Number of stars with determined spectral types. 
'^Extinction values from lBuchholz et al ] l|2009l ). 

Extinction correction using (Afe - 3.0), with Afc from the extinction map in lBuchholz et al ] 1(2001 
"Faintest magnitude bin where the completeness is at least 30%. 
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Table 3. OSIRIS observations of late-type stars 





K' 
IV 




UIIjK^ OHocTi 

\ ) 


■R 
\ ) 


Epoch 




iMa 




"obs 






so- 17 


15.7 


0.02 


-0.09 


0.1 


2004.47 


4 


1.3 


0.17 


11 


20 


c 


SO-6 


14.2 


0.05 


-0.38 


0.39 


2002.69 


4 


2.92 


0.83 


9 


44 


c 


SO-27 


15.7 


0.14 


0.53 


0.54 


2003.34 


1 


2.71 


1.69 


7 


16 


C 


SO-29 


15.6 


0.32 


-0.46 


0.56 


2004.84 


4 


4.63 


1.4 


9 


6 


C 


SO-28 


15.8 


-0.17 


-0.55 


0.57 


2003.53 


4 


1.59 


1.19 


10 


35 


C 


SO- 12 


14.4 


-0.56 


0.39 


0.68 


2002.69 


4 


2.6 


0.51 


4 


33 


C 


SO-13 


13.5 


0.54 


-0.44 


0.69 


2002.69 


4 


2.8 


0.24 


6 


77 


C 


Sl-5 


12.8 


0.34 


-0.92 


0.98 


2002.69 


4 


3.22 


0.64 


4 


45 


c 


Sl-10 


14.9 


-1.13 


-0.03 


1.13 


2002.99 


4 


2.5 


0.86 


5 


25 


c 


Sl-31 


15.7 


-0.98 


0.56 


1.13 


2003.6 


1 


4.62 


1.53 


4 


11 


c 


Sl-6 


15.5 


-0.93 


0.73 


1.18 


2003.34 


1 


3.54 


1.01 


2 


18 


c 


Sl-36 


15.8 


-0.63 


-1.2 


1.35 


2003.58 


1 


2.31 


0.77 


4 


18 


S 


Sl-13 


14.2 


-1.12 


-0.92 


1.45 


2002.91 


4 


2.38 


1.5 


3 


21 


c 


SI- 15 


14.2 


-1.36 


0.52 


1.45 


2002.69 


1 


3.42 


0.27 


4 


38 


c 


Sl-39 


15.4 


-0.55 


-1.37 


1.48 


2003.22 


1 


3.94 


0.66 


7 


16 


S 


Sl-17 


12.5 


0.47 


-1.49 


1.56 


2002.69 


1 


4.33 


0.17 


4 


85 


S 


Sl-44 


15.6 


0.3 


1.61 


1.63 


2003.48 


1 


4.03 


0.4 


5 


22 


N 


Sl-48 


15.5 


-0.62 


-1.53 


1.66 


2004.2 


1 


3.23 


0.87 


10 


32 


S 


Sl-20 


12.7 


0.41 


1.61 


1.66 


2002.69 


1 


4.87 


0.12 


12 


57 


N 


Sl-23 


11.9 


-0.93 


-1.48 


1.74 


2002.69 


1 


3.51 


0.1 


10 


77 


S 


Sl-25 


13.5 


1.65 


-0.63 


1.76 


2002.69 


3 


3.51 




1 


33 


E 


Sl-66 


15.6 


-0.81 


-1.76 


1.94 


2004.87 


1 


1.38 


0.37 


10 


8 


S 


Sl-67 


15.6 


-1.44 


-1.32 


1.95 


2003.56 


1 


2.61 


0.72 


9 


30 


S 


Sl-68 


13.4 


1.86 


-0.61 


1.96 


2002.69 


1 


2.47 


0.52 


4 


50 


E 


S2-12 


15.3 


1.67 


1.17 


2.04 


2003.08 


1 


4.41 


0.16 


5 


18 


NE 


S2-34 


15.4 


1.81 


0.97 


2.05 


2003.43 


1 


3.8 


1 


5 


18 


NE 


S2-11 


12 


1.97 


-0.62 


2.07 


2002.69 


1 


1.33 


0.33 


8 


98 


E 


S2-3 


14.5 


-1.54 


-1.41 


2.09 


2002.84 


1 


2.92 


0.23 


9 


61 


S 


S2-37 


15.8 


0.08 


2.1 


2.11 


2003.91 


1 


3.14 


0.55 


12 


17 


N 


S2-2 


14.1 


-0.55 


2.06 


2.13 


2002.88 


1 


2.26 


0.21 


11 


42 


N 


S2-41 


15.5 


-0.47 


-2.1 


2.16 


2004.41 


1 


2.86 


1.3 


10 


36 


S 


S2-8 


12.2 


-1.84 


0.96 


2.08 


2002.69 


1 


2.82 


0.66 


9 


53 


w 


S2-134 


15.6 


-0.99 


-1.99 


2.23 


2005.79 


1 


3.66 


0.84 


10 


8 


S 


S2-49 


15.5 


-0.83 


-2.1 


2.26 


2004.14 


1 


2.88 


0.46 


10 


19 


S 


S2-47 


14.2 


2.19 


-0.53 


2.25 


2002.75 


1 


2.84 


0.79 


10 


41 


E 


S2-51 


15.9 


0.76 


2.18 


2.31 


2004.88 


1 


4.88 


2.19 


12 


23 


N 


S2-18 


13.3 


-0.98 


-2.14 


2.35 


2002.83 


3 


3.18 


0.17 


10 


57 


S 


S2-23 


14.8 


1.66 


1.77 


2.42 


2002.93 


1 


3.64 


0.38 


4 


36 


NE 


S2-57 


14.7 


-1.21 


-2.12 


2.44 


2003.76 


1 


3.06 


0.3 


10 


45 


S 


S2-26 


13.7 


0.74 


2.42 


2.53 


2002.65 


3 


3.65 


0.34 


12 


43 


N 


S2-62 


15.2 


-1.07 


-2.32 


2.56 


2004.67 


1 


2.89 


0.37 


10 


31 


S 


S2-25 


14.1 


0.74 


-2.45 


2.56 


2003.5 


1 


3.54 


0.21 


8 


91 


SE 


S2-205 


15.8 


0.48 


2.51 


2.56 


2005.54 


1 


3.34 


1.31 


12 


22 


N 


S2-65 


15.8 


2.35 


-1.06 


2.58 


2004.13 


1 


1.96 


0.22 


10 


18 


E 


S2-212 


15.7 


0.18 


-2.59 


2.59 


2006.75 


1 


3.31 


2.01 


7 


12 


S 
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Table 3 — Continued 



Name 


K' 


RA offset 


DEC offset 


R 


Epoch 


Note'' 


Na 






SNR'^ 


Field 






(") 


(") 


(") 






(A) 


(A) 








oZ- 1 U 


14.0 


-Z.U4 


U.41 


9 CKI 
Z.D / 


zUUo.yo 


1 
1 


O. /4 


U.y 


1 1 


90 

zy 


w 


QO Vl 
oZ- / i 


lo.o 


-u.y 


-Z.Oo 


9 flQ 
Z.Do 


oriA/i vo 


1 
1 


1 QO 

i.yy 


U. ID 


1 n 
lU 


oU 


Q 

o 


oZ- 1 Z 


10. i 


-1 .40 


9 9a 

-Z.Zo 


9 71 
Z. / 1 




1 
1 


o.yo 


U.Ol 


1 n 
lU 


Q9 

oZ 


q 

o 


oZ- / o 


10. z 


Z. ID 


-1.D4 


9 79 


ZUUO. 10 


1 

± 


'X HQ 

o.uy 


19 
U. iZ 


1 1 


4y 




oZ- ( o 


1 /I 

14.0 


9 

Z.DO 


-U.oy 


9 7fi 
Z. / O 


OOOO fiO 


1 
1 


9 /lO 

z.4y 


n 7 


1 n 
lU 


oy 


ih 


oZ-oi 


1 Q 1 
lo. 1 


9 fli 
Z.ol 


n 9 
-U.z 


9 Q1 
Z.ol 


OOOO fiO 


1 


z.yo 


n 1 9 
u. iz 


1 n 
lU 


{ 1 


ih 


oZ-oU 


10. D 


9 9Q 
Z. ZO 


1. 1 o 


9 fit; 
Z.oO 


OOOQ ^fl 


1 
1 


9 A 
Z.4 


u.yi 





9n 
zu 


IN IL 


oZ-oi 


10. D 


i.y / 


9 n7 

-Z.U / 


9 Q^J 
Z.oD 




1 
1 


9 Ofi 

Z.yo 


U.Oo 


1 1 
1 1 


99 

zz 




S2-83 


15.6 


2.86 


-0.69 


2.94 


2003.62 


1 


2.17 


1.7 


10 


21 


E 


S2-319 


15.7 


0.99 


-2.79 


2.96 


2005.71 


1 


2.62 


0.81 


8 


14 


SE 


S2-32 


12.4 


1.12 


2.78 


2.99 


2003.05 


1 


3.53 


0.15 


6 


59 


N 


S3-68 


15.8 


-0.35 


3.06 


3.08 


2005.93 


1 


3.96 


0.53 


5 


21 


N 


S3-14 


13.9 


0.1 


3.09 


3.09 


2004.54 


1 


1.98 




1 


18 


N 


S3-16 


15.3 


2.99 


-0.95 


3.13 


2003.76 


1 


5.3 


1.3 


10 


24 


E 


S3-20 


14.9 


1.58 


-2.8 


3.21 


2003.83 


1 


1.64 


0.53 


11 


54 


SE 


S3-6 


12.9 


3.22 


0.01 


3.22 


2002.8 


1 


3.36 


0.26 


10 


77 


E 


S3-109 


15.6 


-3.2 


0.43 


3.23 


2006.75 


1 


4.32 


1.96 


11 


20 


W 


S3-151 


15.7 


1.9 


-2.78 


3.37 


2006.22 


1 


2.99 


1.56 


11 


16 


SE 


S3-32 


15.5 


2.91 


-1.77 


3.41 


2004.51 


1 


3.24 


0.62 


9 


37 


SE 


S3-8 


14 


3.41 


-0.49 


3.44 


2002.81 


1 


2.43 


0.67 


10 


56 


E 


S3- 11 


15.1 


2.96 


-1.93 


3.53 


2004.12 


3 


4.05 


0.84 


4 


35 


SE 


S3-36 


14.7 


3.46 


-0.83 


3.55 


2003.76 


1 


3.65 


0.11 


9 


57 


E 


S3-249 


14.7 


-3.38 


1.32 


3.63 


2006.75 


1 


3.3 


0.34 


11 


33 


W 


S3-37 


15.1 


3.39 


1.35 


3.65 


2004.66 


1 


3.45 


1.23 


5 


17 


NE 


S3-38 


15 


3.69 


-0.11 


3.69 


2005.34 


1 


2.76 


0.92 


9 


33 


E 


S3-39 


13.6 


3.57 


1.1 


3.73 


2002.97 


1 


3.08 


0.53 


3 


32 


NE 


S3-385 


15.3 


3.67 


-1.5 


3.96 


2005.81 


1 


1.32 


0.36 


6 


26 


E 


S4-1 


13.5 


4.02 


-0.37 


4.04 


2003.92 


1 


3.09 


0.25 


9 


102 


E 


S4-2 


13 


3.75 


1.67 


4.1 


2004.01 


1 


3.51 


0.33 


3 


52 


NE 



°'l. Spectral typ e fi rst reported in this w ork. Recent spectral identification: 2. IPaumard et al.l 1I2OO6I ) 3. 
iGenzel et all ||2000| ') 4. ICillessen et~al] ll2009h . 

''Each observation has an integration time of 900 s. 

'^Signal to noise per pixel calculated between 2.212 to 2.218 /im. 
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Table 4. OSIRIS observations of early-type stars 



Name K' RA oflfset DEC oflfset R Epoch Note^ Br 7'' aBr-y ^obs" SNR'l Field 







\ ) 


\ ) 


\ ) 
















SO-2 


14.1 





0.16 


0.16 


2005.61 


2 


5.4 


0.4 


9 


57 


C 


SO-1 


14.7 


0.04 


-0.26 


0.27 


2006.3 


2 


0.5 


0.3 


9 


31 


C 


SO-19 


15.5 


-0.04 


0.35 


0.36 


2006.25 


2 


5.7 


2.3 


11 


16 


c 


SO-3 


14.5 


0.31 


0.12 


0.33 


2006.38 


2 


5.4 


3.6 


9 


21 


c 


SO-5 


15.1 


0.18 


-0.36 


0.4 


2006.32 


2 


7.8 


1.1 


9 


15 


c 


SO-4 


14.5 


0.4 


-0.28 


0.49 


2005.83 


2 


5.4 


6 


10 


1 


c 


SO-7 


15.5 


0.49 


0.1 


0.5 


2006.36 


2 


6.5 


2 


10 


12 


c 


80-11 


15.3 


0.5 


-0.05 


0.5 


2006.73 


2 


10.2 


7.6 


9 


21 


c 


SO-9 


14.4 


0.19 


-0.58 


0.61 


2006.74 


4 


5.4 


2 


7 


38 


c 


SO-31 


15.1 


0.54 


0.44 


0.7 


2006.47 


2 


6.01 


1.54 


6 


21 


c 


SO- 14 


13.7 


-0.77 


-0.27 


0.82 


2005.51 


2 


4.1 


0.7 


7 


76 


c 


80-15 


13.7 


-0.96 


0.24 


0.99 


2005.8 


2 


2.4 


0.5 


7 


67 


c 


Sl-3 


12.1 


0.4 


0.88 


0.96 


2004.98 


2 






5 


79 


c 


Sl-1 


13.1 


1.01 


0.03 


1.01 


2005.2 


4 


1.72 


1 


5 


57 


c 


Sl-2 


14.9 


0.01 


-1.01 


1.01 


2005.93 


2 


4.7 


3.5 


5 


19 


c 


81-4 


12.6 


0.83 


-0.67 


1.07 


2005.84 


1 






5 


52 


c 


81-8 


14.2 


-0.63 


-0.89 


1.09 


2006.17 


2 


5.1 


2.7 


4 


42 


c 


IRS16NW 


10.1 


0.05 


1.22 


1.22 


2004.74 


2 












IRS16C 


9.9 


1.1 


0.51 


1.21 


2003.89 


2 












81-33 


15 


-1.25 


-0.01 


1.25 


2006.26 


1 


5.2 


1.6 


5 


17 


c 


81-34 


13.2 


0.85 


-0.99 


1.31 


2006.38 


1 






3 


52 


c 


81-12 


13.8 


-0.8 


-1.01 


1.29 


2006.21 


2 


3.7 


1.5 


4 


40 


c 


81-14 


12.8 


-1.34 


-0.33 


1.38 


2005.49 


1 


0.5 


0.27 


3 


52 


c 


IRS16SW 


10.1 


1.07 


-0.96 


1.44 


2004.89 


2 












81-19 


13.8 


0.38 


-1.62 


1.66 


2005.45 


1 


7 


1 


4 


37 


s 


81-18 


15 


-0.76 


1.51 


1.69 


2006.56 


1 


1.81 


1 


2 


24 


N 


81-21 


13.3 


-1.66 


0.13 


1.66 


2006.16 


2 






5 


47 


w 


81-22 


12.7 


-1.61 


-0.5 


1.68 


2005.65 


2 


-1.61 


0.49 


3 


61 


c 


81-24 


11.6 


0.71 


-1.62 


1.77 


2005.64 


2 


2.84 




1 


74 


SE 


82-5 


13.3 


1.91 


-0.81 


2.07 


2005.68 


1 






7 


58 


E 


IRS16CC 


10.6 


1.99 


0.56 


2.07 


2004.57 


2 












82-4 


12.3 


1.47 


-1.47 


2.08 


2006.04 


2 


1.94 


0.19 


11 


149 


SE 


82-6 


12 


1.61 


-1.35 


2.1 


2004.69 


2 


2.04 


0.16 


7 


103 


SE 


82-7 


14.1 


0.96 


1.84 


2.08 


2006.15 


2 







12 


54 


N 


IRS16SW-E 


11 


1.86 


-1.14 


2.18 


2005.52 


2 












IRS33N 


11.4 


-0.05 


-2.21 


2.21 


2006.19 


2 


3.1 


0.4 


10 


104 


s 


82-17 


10.9 


1.29 


-1.88 


2.28 


2005.14 


2 


3.19 


0.07 


11 


125 


SE 


82-50 


15.5 


1.68 


-1.52 


2.26 


2005.38 


1 


8 


3 


11 


16 


SE 


82-22 


12.9 


2.3 


-0.24 


2.31 


2005.33 


1 


-0.73 


0.5 


10 


78 


E 


82-19 


12.6 


0.42 


2.31 


2.34 


2005.53 


2 


2.04 


0.5 


12 


68 


N 


82-21 


13.5 


-1.66 


-1.64 


2.34 


2006.34 


1 


1.2 


0.7 


8 


71 


s 


82-55 


15.4 


0.88 


-2.19 


2.36 


2006.53 


1 


3.2 


1.2 


11 


29 


SE 


82-58 


14.1 


2.14 


-1.16 


2.43 


2005.74 


1 


4.7 


3 


7 


41 


E 


82-59 


15.4 


0.79 


-2.34 


2.46 


2006.48 


1 


5.6 


2 


10 


25 


SE 


82-61 


15.5 


2.37 


-0.67 


2.47 


2006.56 


1 


6 


3 


10 


16 


E 
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Table 4 — Continued 



Name 


K' 


RA offset 
(") 


DEC offset 
(") 


R 

(") 


Epoch 


Note'' 


Br 7*^ 

(A) 


(A) 




SNR<1 


Field 


S2-74 


13.3 


0.15 


2.78 


2.78 


2006.16 


2 







12 


81 


N 


S2-77 


13.8 


-1.78 


-2.19 


2.82 


2006.16 


1 


1.4 


0.5 


9 


54 


S 


S2-76 


15.2 


-0.23 


2.81 


2.82 


2006.76 


1 


1.18 


1 


12 


38 


N 


S2-317 


15.6 


-0.74 


-2.86 


2.95 


2006.4 


1 


5.1 


1 


3 


21 


S 


IRS16NE 


9.1 


2.87 


1.01 


3.04 


2003.62 


2 












S3-2 


12.1 


3.06 


0.53 


3.11 


2005.43 


1 


0.49 


0.04 


3 


47 


NE 


S3-3 


15.1 


3.07 


-0.67 


3.14 


2006.47 


1 


1.9 


2 


10 


21 


E 


S3-5 


12.2 


2.94 


-1.18 


3.16 


2005.66 


2 






10 


87 


E 


S3-7 


13.9 


1.91 


-2.61 


3.24 


2006.34 


1 







11 


76 


SE 


S3-30 


12.8 


1.64 


-2.96 


3.38 


2006.39 


2 


1.63 


0.44 


8 


56 


SE 


S3-34 


13.5 


3.22 


-1.26 


3.46 


2006.29 


1 


2.12 


1.28 


9 


46 


E 


S3-190 


14.1 


-3.18 


1.45 


3.49 


2006.47 


1 


1.2 


1.1 


11 


55 


W 


S3-10 


12.4 


3.33 


-1.14 


3.52 


2005.99 


2 






10 


81 


E 


S3-12 


12 


2.37 


-2.73 


3.61 


2001.39 


1 







11 


114 


SE 


S3-314 


15.4 


3.81 


-0.12 


3.82 


2006.43 


1 


6.8 


2.6 


9 


26 


E 



'^l. Spect r al typ e first reported in this work. Recent spectral identification: 2. IPaumard et al.l ||2006| ) 3. 
iGenzel et all bOOfj) 4. ICillessen et~al] boo^ ). 

'^Each observation has an integration time of 900 s. 

Stars with no reported equivalent widths have no detectable Br 7 lines (see Sectional, with the exception of the 
IRS 16 sources, which all have strong Br 7 emission lines from strong winds. 

"^Signal to noise per pixel calculated between 2.212 to 2.218 /^m. 



Table 5. OSIRIS spectra with unknown spectral type 



Name 


K 


RA ofTset 

{ ) 


DEC offset 

( ) 


R 

( ) 


Epoch 


SNR' 


SO-35 


15.3 


0.06 


0.97 


0.97 


2003 77 


16 


Sl-26 


15.6 


-0.92 


0.39 


1 


2003 07 


12 


Sl-29 


15.4 


1.08 


0.17 


1.09 


2005 09 


11 


Sl-32 


15.3 


-0.96 


-0.64 


1.15 


2003.21 


20 


Sl-7 


15.9 


-1.03 


-0.53 


1.16 


2003 86 


13 


Sl-85 


15.5 


0.92 


-0.82 


1.23 


2006 75 


4 


Sl-40 


15.9 


-1.39 


-0.62 


1.52 


2004 23 


7 


Sl-41 


15.8 


0.99 


1.15 


1.52 


2003 69 


4 


Sl-42 


15.9 


0.94 


1.26 


1.57 


2004.04 


9 


Sl-45 


15.3 


-1.27 


1.04 


1.65 


2003 04 


13 


31-47 


15.6 


-1.58 


0.46 


1.65 


2003.14 


1 


Sl-50 


15.4 


1.5 


0.68 


1.65 


2003.7 


9 


Sl-128 


15.8 


1.33 


-0.99 


1.66 


2006.55 





Sl-51 


15 


-1.65 


-0.18 


1.66 


2002.99 


10 


Sl-52 


15.2 





1.67 


1.67 


2003.15 


28 


Sl-53 


15.3 


1.67 


-0.12 


1.67 


2003.12 


11 


SI- 54 


15.5 


-1.52 


0.74 


1.69 


2004.25 


3 


Sl-55 


15.5 


1.59 


0.6 


1.69 


2004.17 


7 


Sl-56 


15.8 


-1.11 


1.35 


1.75 


2003 85 


11 


Sl-61 


15.9 


-1.49 


-1 


1.79 


2003.96 


6 


Sl-62 


15.3 


0.5 


1.75 


1.82 


2003.21 


19 


Sl-59 


15.5 


0.02 


1.82 


1.82 


2003 33 


27 


SI- 147 


15.9 


1.18 


-1.47 


1.89 


2006 75 


7 


81-64 


15.5 


0.64 


1.81 


1.92 


2003 08 


13 


Sl-159 


15.9 


1.3 


-1.43 


1.93 


2005 94 


9 


Sl-65 


15.7 


1.46 


1.27 


1.93 


2003 85 


5 


S2-89 


15.9 


1.07 


-1.7 


2.01 


2006 75 


5 


S2-42 


15.6 


0.46 


2.1 


2.15 


2004 02 


23 








n 




ZUUO. ZrO 


12 


S2-63 


15.4 


-0.65 


2.48 


2.56 


2003.98 


12 


S2-219 


15.8 


-1.6 


-2.07 


2.62 


2005.34 


8 


S2-261 


15.6 


-2.61 


1 


2.8 


2006.04 


11 


S2-268 


15.5 


-2.78 


0.49 


2.82 


2006.32 


13 


82-304 


15.9 


2.52 


1.47 


2.91 


2006.39 


12 


S2-306 


15.6 


-0.48 


-2.89 


2.93 


2006.75 


17 


S2-84 


15.6 


1.65 


-2.47 


2.97 


2004.29 


26 


S3-48 


15.9 


-3 


0.37 


3.02 


2006.32 


15 


S3-50 


15.7 


-1.87 


-2.38 


3.03 


2006.75 


13 


S3-51 


15.1 


-0.15 


3.02 


3.03 


2005.88 


15 


S3-62 


15.7 


2.68 


-1.48 


3.06 


2006.4 


28 


83-65 


15.7 


-1.27 


-2.8 


3.07 


2005.72 


11 


S3-4 


14.7 


3.08 


-0.52 


3.13 


2002.91 


28 


S3-90 


15.8 


2.74 


-1.56 


3.15 


2005.21 


30 


S3-91 


15.7 


-1.73 


-2.67 


3.18 


2006.75 


8 


S3-116 


15.9 


2.25 


2.37 


3.27 


2006.11 


9 



-35- 



Table 5 — Continued 



Name 


K' 


RA offset 
(") 


DEC oflfset 
(") 


R 

(") 


Epoch 


SNR* 


S3- 24 


15.6 


3.26 


0.41 


3.29 


2004.87 


4 


S3- 167 


15.5 


-3.08 


1.42 


3.4 


2006.75 


16 


S3- 169 


15.5 


-3.39 


0.38 


3.41 


2006.75 


27 


S3-31 


15.5 


3.4 


0.38 


3.42 


2004.42 


7 


83-33 


15.4 


3.32 


-0.85 


3.43 


2004.25 


16 


83-286 


15.3 


3.42 


-1.54 


3.75 


2004.66 


13 


83-319 


15.4 


3.55 


-1.42 


3.83 


2006.39 


17 


S4-47 


15.9 


2.56 


-3.23 


4.12 


2006.75 


11 


S4-67 


15.7 


2.8 


-3.05 


4.14 


2006.75 


7 



^Signal to noise per pixel calculated between 2.212 to 2.218 /im. 



